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CN ■ rotating strings and D5 brane, which is introduced as the baryon vertex. Our model is 
! based on the type IIB superstring theory with the background of asymptotic AdS^ x 5"^. 
This background is dual to a confining supersymmetric Yang-Mills theory (SYM) with 
gauge condensate, (F"^), which determines the tension of the linear potential between 
the quark and anti-quark. Then the slope of the meson trajectory (a^) given by 
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^ [ this condensate as a^,/ = 1/^'k{F'^) at large spin J. This relation is compatible with 
\ the other theoretical results and experiments. For the baryon, we show the importance 
of spinning baryon vertex to obtain a Regge slope compatible with the one of and 
A series. In both cases, mesons and baryons, the trajectories are shifted to large mass 
side with the same slope for increasing current quark mass. 
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1 Introduction 



Among many observations of hadronic phenomena, the Regge behavior, J = ao + 
a'^^M"^ , for the hadronic mass M and its spin J has imphed the existence of strings 
which connect constituent quarks. The string theory, however, has been estabhshed as 
the quantum theory of the gravity in ten dimensional (lOD) space-time. On the other 
hand, SU{'i) Yang-Mills theory has been established as the basic one to form hadrons 
in 4D space-time. So the two theories seem to be separated as different one. 

After being uncovered these facts, new paradigm in the research of the strong cou- 
pling Yang-Mills theory has been opened by the conjecture known as gauge/gravity 
duality or AdS/CFT correspondence [H Ej |3]. Many trials to understand strongly 
interacting gauge theories have been revived within this paradigm. The approach is 
based on the lOD string theory, which describes gravitational theory. Using the low 
energy effective action of this theory, many classical solutions have been studied as the 
dual of the 4D Yang-Mills theories. This approach is also called as the holography 
since the 4D gauge theory is reflected in the higher dimensional gravity. 

Here we study the Regge behavior of hadrons from the gauge/gravity duality. Two 
possible ways are considered in this approach. One way is to examine the fluctuation 
modes of probe brane embedded in the supergravity background dual to the confining 
theory. The approaches in this direction have been performed from five dimensional 
effective gravity [H El E]. The string theory is however not yet solved in curved space- 
time, then it is still an open problem to introduce higher spin states in a strict sense. 

Another way is to examine classical string-configurations in the supergravity back- 
ground dual to the considering Yang-Mills theory. In this case, higher spin meson-states 
are introduced through the rotation of the strings, whose two end points are on the 
flavor brane introduced as a probe. While such states are not quantized, we can see 
the relation of mass and spin for mesons in terms of these solutions at least for large 
spin. Here we study the Regge behavior through classical solutions as performed in 
the references [3 El El EDI [IH [El [13] . 

In the M=2 supersymmetric case with flavor brane(s) [71 [TTl [12], the background is 
given by AdS^ and Regge-like behavior is observed at small J as 



where and A denote the current quark mass and the 'tHooft coupling. In this case, 
however, the Regge behavior is not observed at large J because the potential between 
quark and anti-quark is not linear at large distance. The behavior in the region of 
small J is common to the theory dual to the background which approaches to AdS^ in 
the UV limit or near the boundary. Then this behavior corresponds to the behavior of 
the CFT. Here skipping this CFT behavior at small J, we concentrate on the Regge 




(1.1) 
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trajectory given at large J in the confining theory. 

In the background dual to the confining theory, we expect a linear Regge trajectory 
at large J. In the D4 branes background with fiavor D6 probe brane model, such 
behavior has been observed P, [101 [H] at the limit of rriq — > 0, 

.^„'(|lfM% ,1.2) 

where a', Rda and Uh represent the string tension, the radius of of bulk space-time 
and the scale of compact S*^ of world volume of D4 brane respectively. It should 
be noticed that main parameters of the bulk lOD background are remained in this 
formula. For finite m^, however, the trajectory deviates from the linear behavior. This 
deformed behavior has been discussed related to the heavy quarkonium [10]. In any 
case, as a common property of this kind of classical model, J arrives at zero for M = 0, 
and then the intercept of the trajectory is always zero. In order to have a non-zero 
intercept, quantum corrections would be necessary as pointed out in [6l |8]. Here we do 
not extend our analysis to such quantum version. We study the Regge trajectory by 
using only the classical solutions. 



In our approach, the lOD background configuration, which is dual to a confining su- 
persymmetric Yang-Mills theory (SYM), is set as the solution of type IIB supergravity 
with five form field fiux, dilaton and axion [T^ [T5] . In this model, the tension tm of 



the linear potential between the quark and anti-quark is given as tm = y{F^)/'2 by 
the gauge condensate (F^) ^7\. The mesons and baryons are constructed by embed- 
ding probe D7 brane (fiavor brane), fundamental strings (F-strings) as quarks and D5 
brane as baryon vertex. The D7 brane is introduced to incorporate fiavor quarks. The 
F-strings could end on this D7 brane. These are all introduced as probes, so they do 
not change the background configuration. That is, the vacuum of the gauge theory is 
not altered by these branes. 

The meson with higher spin is given by a rotating F-string whose both ends on the 
D7 brane. Its stable configuration is obtained by solving the equation of motion for the 
Nambu-Goto action of the F-string in the bulk background. In this case, appropriate 
boundary conditions for strings on the D7 brane are imposed. Since the mass and the 
spin of the meson are not quantized as mentioned above, then they are continuous. At 
large J, as expected, we find a linear trajectory, 

J= / , (1.3) 



with the slope parameter = l/{2nTM) = 1/tx\J{F^. We should notice that the 
above formula is expressed only by the quantity of the gauge theory side without any 
parameter of the bulk side. This result is consistent with the picture obtained from 
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naive quark model. That is, the rotating quark and anti-quark are connected by a 
string with the tension of tm- On the other hand, at small J, we could find the 
behavior found in the case of the AdS^ background mentioned above. The behavior in 
this region depends on the current quark mass ruq, which is determined by the profile 
function of the embedded D7 braneo As mentioned above, this meson configuration 
gives J = at M = 0, then the intercept is zero. 

In the next, we extend this approach to the baryon trajectories, where the baryons 
are constructed according to the model proposed before [l8]-[23]. To make baryons, 
we must introduce D5 brane as the vertex which attracts Nc F-strings. The F-strings 
connect the D5 and D7 branes, and then both F-strings and D5 brane contribute 
to the spin and mass in the case of the baryon. While the D5 vertex could have 
various configurations, for the simplicity, we examine here two typical configurations 
given in [23]. One is called as "point vertex", which does not extend in the real three 
dimensional space. So it looks like a point in our space. Another one is called as "split 
vertex" , which extends in one direction in our three space like a string. From the two 
end point of this vertex, F-strings are going out toward the D7 brane. 

In the case of the point vertex, spin comes from the rotating F-strings only since 
the point vertex couldn't generate rotation. As a result, we find that the Regge slope 
of this baryon is given by a'^ = 2a^^/(A^c(l -|- /?)) with a positive constant (3. In this 
case. Then it may be difficult to realize the desirable slope parameter in this case for 
iVc>3. 

On the other hand, for the case of split vertex, the smaller the energy of the state 
is, the shorter the attached strings become, then the configuration may be dominated 
by the vertex. The vertex volume doesn't shrink to zero and remains with a finite size 
then with a finite energy. This means that the lowest energy state is finite contrary to 
the mesons. Another point to be noticed is that the split vertex rotates to generate 
the spin. So we estimate the spin and the mass of the baryon by the vertex only in 
this case. Then the Regge slope of this configuration is given by the tension of the 
split vertex. We find a reasonable slope parameter for the case of A^^^ = 3 and 4. It is 
favorite to consider the case of A^ = 3, but a problem to be solved is remained in the 
case of A^c = 3 as discussed below. Here we add the analysis for A^c = 4 since this case 
is free from such a problem. 

The outline of this paper is as follows. In the next section, the bulk solution for our 
holographic model is given. In the section 3, meson trajectories are studied by using 
Nambu-Goto action. Introducing D5 brane, the baryon trajectories are studied in the 
section 4. The summary and discussions are given in the final section. 

^This is given by the value of the profile function at UV limit of the fifth coordinate (r). However, 
in the present supersymmetric case, the profile function is a constant. 
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2 Holographic Model 



Here we give the ingredient of our holographic model for supersymmetric Yang-Mills 
theory with quarks to form meson and baryons. 

2.1 Bulk Background 

We consider lOD IIB model retaining the dilaton $, axion x ^i-nd self-dual five form 
field strength F^^y Under the Freund-Rubin ansatz for F(5), F^^...^^ = — a/A/2 e^j...^g , 
and for the lOd metric as M5 x 5"^ or ds^ = g^^dx^dx'^ + gijdx^dx^ , the solution given 
below has been found [HI [15]. Where (/x, z/) = ~ 4 and (z,j) = 5 ~ 9. The five 
dimensional part (M5) of the solution is obtained by solving the following reduced 5D 
action, 

3=^2/ d'^^ {r + 3A- + le'^idxf) , (2.1) 

which is written in the string frame. 

The solution is obtained under the ansatz, 

X = -e^" + Xo , (2.2) 
which is necessary to obtain supersymmetric solutions. The solution is expressed as 
ds'^Q = G ]\^jj\[dX^^ dX^ 

Then, supersymmetric solution is obtained as 

e* = l + ^, A = l, (2.4) 

where M, iV = ~ 9 and = 71/2 = (AirgsNca'^/^ = (Aa'^)^/^ and A = ATrg^Nc 
denotes the 'tHooft coupling. The dilaton is set as = 1 at r — )■ 00, and the parameter 
q corresponds to the vacuum expectation value (VEV) of gauge fields strength [15] of 
the dual theory. Then this solution is dual to the four dimensional A/'=4 SYM theory 
with a constant gauge condensate. Due to this condensate, the supersymmetry is 
reduced to M=2 and then the conformal invariance is lost since the dilaton is non- 
trivial as given above. As a result, the theory is in the quark confinement phase 
since we find a linear rising potential between quark and anti-quark with the tension 
^/(27ra'i?2) [lll[T5l[I7]. 

Furthermore, we can see that the space-time is regular at any point. In the ultraviolet 
limit, r — J- 00, the dilaton part e'^ approaches to one and the metric (12. 3p is reduced to 
AdS^ X . On the other hand, the dilaton part e*^ diverges in the infrared limit r — j- 0, 



e*/2 [^A\r) [-de + {dx'f) + ^dr^ + R^d^A . (2.3) 
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so that one may expect a singularity at r = 0. However there is no such a singular 
behavior. This is assured by rewriting the metric (12. 3p in terms of new coordinate 2;, 
where z = R'^/r. Then we obtain 




In the infrared limit 2; — )■ 00, we have 

Therefore we find lOD flat space time in this limit and no singular point [Til [T5] . 



2.2 D7 brane embedding 

In this background a probe D7 brane is embedded to introduce Af=2 hypermultiplet, 
which corresponds to a string connecting the stacked D3 branes and the probe D7 
brane. This is regarded as the quark hereafter. The world volume of the D7 brane is 
set by rewriting the extra six dimensional part of (12. 3 p as 

^dr^ + RW, = ^{dr^ + r^dnl)=^(±dX^^' 




dr]^ + 7]''dnl + Y.dX' ] . (2.7) 

i=8 ) 

Then the world volume coordinate of D7 brane is taken as the four dimensional space- 
time and extra four dimensions, (77, ^3). Then the induced metric on the D7 brane 
is expressed as 

R^ 



dsi 



"''R^ I ^^A\t) (-de + {dx^) + 4 ((1 + ^\v?)dr,' + v'dnt) \ . {21 



where = 77^ + w{riY and w'{ri) = dnw{ri) Here the D7 brane is embedded under the 
ansatz, 

{X^f + {X'^f = w^{r]). (2.9) 

The solution of ^(77) is obtained as w =constant in the background considered here 
[17]. While, for the simplicity, the solution is given as (X^,X^) = (w,0) in [17], we 
can choose it at any point of the circle [X^Y + {X'^Y = ^ since the background is 
symmetric in the (X^,X^) plane. Within the probe approximation, A'^^ << N^., Nf 
number of D7 branes may be put on this circle. We need here at least two D7 branes, 
Nf = 2, to obtain the split vertex given in the section 4. Anyway, for this solution, 
we can make a simple image for the baryon configuration as shown in the section 4. 
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Remembering the relation, z = E? /r, the distance between D7 and D3 branes is given 
by a constant rm = R'^/zdj, then the current quark mass is obtained as 



rD7 



m„ = -^ = ^-— . (2.10) 



2.3 Rotating string 



Mesons are given by a open string which two end points are on the D7 brane. The 
baryons are made by Nc strings, whose one end is on the D7 brane and the other end 
is on a vertex for every string. The vertex is given by a probe D5 brane which wraps 
on of extra dimension[T8] in this modeL Here the rotating string is formulated as 
follows. 

In the following analysis we change the coordinate r as R^/r = z, then we rewrite 
the metric (12. 3p as follows 



dsj, = e"/^^ {A'iz) (-de + {dx'f) + dz" + z^dn''^ . 



(2.11) 



Further, the metric for the string which rotates around the axis is rewritten by 
cylindrical polar coordinates as. 



dsl^ = e"/^^ [A'iz) (-de + dp" + pHe^ + rfx^) + dz' 



(2.12) 



Taking the string world sheet as (r, a) = {t, z) and the ansatz, p = p{z) and 6 = out, 
the induced metric is given as 



grr = e''/'^A'{z) (-l+cV 



9o 



,<E>/2 



1 



where prime denotes the derivative with respect to z. Then we have 

^string = JdtC=-^J dtdze''/'A\z)^^{l - a;V) {p" + 
From this, the spin Js and the energy Es of this string are given as 

i?2 



(2.13) 



(2.14) 



J. 



dC 



1 



dC . 
u- £ 

OOJ 



- dze^/^A\z)-^up' 



p" + A-\z) 



uj'^p^ 



2na' 



dze^/'A'iz)^. 



p'' + A-'^{z) 



\ 1 



uj'^p^ 



(2.15) 



(2.16) 
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These are estimated by giving appropriate solutions for the corresponding strings in 
mesons and baryons as shown below. 

Here we emphasize that the D7 and fundamental strings introduced here are all 
treated as probes, so they do not alter the bulk background and their embedded con- 
figurations. We only impose consistent boundary conditions at their connecting points 
in constructing the hadron configurations. As for the D5 brane, we explain it in the 
section 4, where the baryons are discussed. 

3 Regge trajectory of Meson 

The meson with higher spin is constructed by rotating F-string whose end points are 
on the D7 brane. The configuration of such a meson state is obtained by solving the 
equation of motion for the Nambu-Goto action fl2.14p with an appropriate boundary 
condition at the D7 brane. Using this solution, we find the relation of the spin fl2.15p 
and the mass ( 12.16^ to obtain the Regge trajectory. 

3.1 Rough Estimation of Slope parameter for small niq 

It is difficult to obtain analytic solution for the string. So before giving the numerical 
estimations, we give a rough estimation of the relation J (x E"^ in the case of light quark 
mass. At the small quark mass limit, rriq — )■ 0, the speed of the quark approaches to 
the light velocity, upb — )■ 1, where pb denotes the value of p at the end points of the 
string. They are on the D7 brane. In this limit, the integrations of fl2.15p and fl2.16p 
are approximated as follows. 

At first, we rewrite f l2.15p and fl2.16p by changing the integration variable from z to 
p as 



Here we can regard M{p) and I as the mass density and the inertial moment of the 
rotating string respectively. Since M{p) depends on p through z{p), which is obtained 
as the solution of the equations of motion, it is hard to perform the integration in (13. ip . 
In the present case, we can however estimate it by a reasonable approximation, which 
is ensured by the numerical estimations given later. 
M(p) is written as 




(3.1) 




(3.2) 




(3.3) 



(3.4) 
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We notice that Mgta is the integrand in estimating the potential of the static quark and 
the anti-quark. As we have previously shown, we find a linear potential in our present 
model at large distance between the quark and the anti-quark. This implies, for large 
Pb, that we obtain the following result 

^sta = 2 / rfpMsta = 2rA/P6 , (3.5) 

Jo 

where the tension, tm, is given in our model as follows [1 



This result is understood as follows. For large p^, the string shape approaches to 
the U shape with a long bottom horizontal-line, which is the main part of the string 
configuration. In the range of this part, / z = r is small and almost constant. Then 
in this region, we can set as 

for the supersymmetric case (12.41) . which is written also as, 

e* = l + gV, g = -|^, A = l. (3.8) 

Hereafter, we consider this solution for the simplicity. In this approximation, the both 
side parts of U-shaped string configuration are neglected. They corresponds to the 
effective quark and anti-quark masses. So the approximation is applied to the case of 
small quark mass. 

Here we notice the following relation 



g = vr2(F2)Aa'\ (3.9) 

where A = AnggNc and (F^) denote the 'tHooft coupling and the VEV of the gauge 
field condensate. Then, remembering = Aa'^, we find 

1 



TM = -^{F^) (3.10) 

This implies that the tension does not depend on the other parameters, R and a', and 
then it is determined only by the mass scale, (-F^), of the confining gauge theory. We 
notice here another mass scale, i?, which determines the meson mass obtained from 
the D7 brane fiuctuations [16]. We also observe the Regge behavior with this mass 
scale when we restrict the region to the small mass and to the finite quark mass case. 
In this case, the trajectory is given by (11. ip as mentioned in the Introduction. This 
behavior is seen as a common properties of the meson spectra through the probe brane 
fiuctuations. 
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Then A^sta can be approximated by the QCD-string tension tm for large J and small 
quark mass. Under this approximation, Em and Jm are estimated as follows, 



rPb / 1 fpb 0^ 

Em = 2tm / dpJ- — , / = 2rM / dp ./ . , • (3.11) 

Jo V 1 ~ ^ P "'0 V 1 — ^ P 

and using the approximation up^ ^ 1, we get 

TT 

Em = T^TMPb , Jm = -^mPb' ■ (3.12) 



As a result, we obtain 



Jm = tt—Em^ ■ (3-13) 
This implies that the slope of the meson trajectory is given by 

The result f l3.13p is interesting since its form is parallel to the one obtained in the flat 
bulk-space analysis [8], where the Regge slope is obtained as 

We flnd that our result is obtained by replacing the F-string tension Tq by the QCD- 
string tension tm in the above formula. This point is well understood from the 
fact that we flnd lOD Minkowski space-time in the IR limit of our model as shown 
above. However, we should notice that the scale of the coordinates are redeflned as 
{q^/^/R)x>' x^. 



3.2 Quark mass effect 

In the next, we consider the case with a flnite m^, which can not be negligible. We 
rewrite the lagrangian f l2.14p according to the formulation performed in [10] as follows, 

£= L{z)dz= r' L{z)dz+ r L{z)dz + H L{z)dz , (3.16) 

J Z— J Z— J Zi J Zr 

where 

L{z) = -^P(^)^(l-u;V) (P" + 1), (3.17) 
and ^ 

P{z) = e''^^A\z)^. (3.18) 
z^ 

The coordinates z± denotes the two end points of the F-string at the D7 flavor brane. 
The region (i) zi < z < Zr represents the bottom range of the U-shaped string con- 
flguration, which is considered in the above subsection. Then, the two parts, (ii) 



9 



< z < zi and (iii) Zj. < z < corresponding the rotating quark and anti-quark 

respectively. Since the configuration is symmetric, in the present case, then they are 
equal. And they are neglected in the above as very small quantity. Here we estimate 
them. 

We solve the string equations in terms of the above lagrangian, ( I3.16p . The solution 
has the U-shape, so we approximate as 



|p'(^)|»l, 

for region (i) 

for region (ii) and (iii). Then the variation with respect to p is given as 



(3.19) 



(3.20) 



dz{-d,\Piz)Jl 



9 9 



P' 



V2 



P(z] 



P 



/2 



VI - uj-'p' 



5p 



2'Ka' 



2Tia' 



dz P(z) 



P(z] 



-oo'^p 



uj'^p^ 



V2 



5p 



2+ 



dz P(z) 



-oo'^p 



-,5p 



(3.21) 



Vl — W^P^ Jzr \/l — Uj'^P^ 

The second boundary terms come from the partial integration of the variation for the 
first term of fl3.16p . This must be cancelled out with the terms in the third parenthesis 
of fl3.2ip since the first term vanishes due to the equation of motion. Noticing the 
symmetric configuration, we find the next condition. 



1 2 2 



Zna — — rm, 
Piz*] 



(3.22) 



2na' 



dz P(z) 



(3.23) 



where rUg denotes quark mass. Here, further we set as zi = Zr = z* and p{z*) = pb- 
The negative sign of the above fl3.23p comes from p'{zj) ~ — oo. We notice that the 
position of the D7 brane is given by rj = / z±. Since z^ < Zr, we have positive niq. 
The same relation between ph and nig is obtained in from the other side condition. 
This relation is consistent with the fact that the speed of the massless quark is one, 
namely uph = I for nig = 0. 

Then, under the condition of small nig but non-zero, we estimate Jm and Em ac- 
cording to (12.151) and (I2.16P by separating the region of z as above. Using the condition 
given above at the separating point, we obtain 



E 



M 



7rrjv/(pfe 



l/2^ 



(3.24) 
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JM = \'nrM[pl + 2^pl'') , (3.25) 



where 



2 2Txa'm„ 



Then Jm is expressed by Em by ehminating p^. It is further expanded as a series of 7 

as 

The second term is negative and proportional to m^. Then we would find a shift of the 
trajectory to the right in the J — E"^ plane as the quark mass effect. We should notice, 
however, that the slope given at large Em is not changed by rriq. 

3.3 Numerical estimations 

Equations of motion for SUSY background 

The values of these quantities are obtained by introducing the D7 brane to consider 
quarks with the finite mass. Here, we consider the supersymmetric case, (12. 4p or (13. 8p . 
In order to solve the string equation, it is convenient to use the reparametrization in- 
variant formalism. The solutions p{z) are multi-valued for z, however the configuration 
of one solution is given by a continuous curve. So the solution can be expressed by one 
parameter, s as given here. The Lagrangian is written by using the parameter s as 

C = I dsL = £ dse'/'A\z)^^{l - uj^p^) (p^ + z^A-\z)) . (3.28) 

where dot denotes the derivative with respect to s. si and Sr are defined as 

z{si) = z{sr) = zdt , p{si) = p{sr) = . (3.29) 

The first two equations means that the end points of the string are on the D7 brane, 
then the second equations are coming from the Dirichlet condition p'{ziji) = for the 
string end points on the D7 brane. We impose these boundary conditions in solving 
z{s) and p(s) by using (I3.28p . 

The equations to be solved are obtained as follows. Introducing the canonical mo- 
mentum as, 

dL dL , , 

Pp=gr, P.= gr, (3.30) 



we have the Hamiltonian 



H F 
H = 2-, A= (3.31) 

^ ,/p2 + PA-^{z) 
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F = e''/^A\z)^^l-u'^p\ (3.32) 



z 

H = 1{pI + pIA\z)-F') . (3.33) 
Then the Hamilton equations are obtained from H instead of H for the simphcity, 

p = p,, z=p,A\z), (3.34) 

Pp = -ujpQ{z), p, = -p^A[z)^— + -[l-uj p j , (3.35) 

and 

Q = e'A\z)-^. (3.36) 

In solving the above simultaneous equations (13.34^ and ( 13.35^ . the boundary condi- 
tions fl3.29p are imposed as mentioned above. And we need the values of the parameters 
used in the equations to get numerical results. 

About the value of A = R^/a'"^ 

We expect holographic approaches to be useful at large A = R^joi'^ and large Nc- 
Although, in our real world, iVc = 3 and this is not large enough, we dare to compare 
our results with the experiments obtained up to now. We know the p meson trajectory, 
which is approximated by J = 0.53 + O.SSi?^. So, from (13.141) . we can set as 

^Regge ^ 1 ^ Q gg (GeV-2) . (3.37) 

This implies 

= 0.36 (GeV^) . (3.38) 

This is our prediction for the gauge condensate. As for the value of this quantity, 
many attempts to estimate it have been done in 4D gauge theory side. We compare 
our result with the one given in [2i^ 



ig^uNcp2^ = 0.14 (GeV^) , (3.39) 

according to the lattice simulation. Here the value is obtained for quenched approxi- 
mation and for = 3. Our calculation is also quenched for the quarks, then they can 
be compared. Noticing g\^Nc = AirggNc = A, we find 

A = 42 . (3.40) 

Then A is not small, thus it may be reasonable to compare our result with the quantities 
in our real world. 
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In applying the above result (13.401) . we should notice that our model is based on the 
self-dual gauge field configurations which satisfies FF = FF. Then we should study 
also the value of FF to check the validity of using the result given in the lattice theory, 
where the supersymmetry is however lost. As in the case of FF, we have searched for 
simulations of FF in lattice theories, but we could not find any definite result for this 
quantity. The lattice formulation might be unsuitable to study topological quantities 
like FF. 

Even if the self-dual relation for the gauge fields in the background were largely 
broken in the real QCD, it would be still meaningful to show our analysis since we 
expect that our main results would survive as being approximately true in the real 
QCD. Although this point is open, our results are not so sensitive to the value of A 
given above as shown below. 

Numerical results 




Fig. 1: The numerical result of J - E'^ relations for q = l/(.88^ x 42) = 0.0306 and R^/a' = 
6.48. The curves (a) and (b) denote the results for I/z-qj = 0.005, and 0.50 respectively. 
Here the curves are almost the same for 20 < A < 42. And the small (large) dots denote the 
p (tt) meson series. 



In obtaining the numerical values of J and E, we must fix the position of the D7 
brane or the quark mass X^^"^ / {2'kzd'j) = m^, and other parameters, g ( or g) and / a' 
(or A). The latter two are fixed by neglecting the quark mass and using our formula, 
experimental p meson trajectory and a result of lattice simulation. The results are 
given in the previous subsection (3.2). 
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In FigUl two numerical results are shown for q = 0.0306 and R^/a' = 6.48 with 
the mass spectrum of p and tt meson series. The two curves, (a) and (b) are for 
l/^^D? = 0.05 and 0.5, respectively, where 1/2:07 denotes the position of the D7 brane 
and it is equivalent to the quark mass m^. For very small ruq the trajectory is linear 
in almost all region and is consistent with fl3.14p . Further we find the expected shift 
to the right of the trajectory for the large rUq case, I/^d? = 0.5. As for the quark mass 
dependence, we can see that the heavier the quark mass becomes, the larger the energy 
is for the same J. This is a reasonable result. The value of A is set as 42 in the figure, 
we notice however that the results are the same for wide range of A, 20 < A < 42. 

We notice one more important point that the slopes of the two curves are different 
at small mass region. While, for any quark mass, we could find tha same slope given 
by f l3.14p at large meson mass as discussed above, we can see the difference of the 
slope, which depends on the quark mass m^, in the region of small meson mass. This 
behavior in small meson mass region is consistent with the trajectory given by (11. ip . 
In this sense, then we could see two regions where the mass scale of the string tension 
is different. 

However we should notice that both of the curves approach to the origin for J — )■ 0. 
This implies that the intercept of the trajectory gives zero in the present case. This 
is understood as follows. In our calculation, the classical solutions of the Nambu-Goto 
action are used, and they are U-shaped and give positive E. When E vanishes, the 
string configuration also disappears since the calculation is classical. Then, J vanishes 
simultaneously for E = 0. This implies that the intercept should be zero in this 
calculation. One way to get finite intercept, in this approach, is to add the quantum 
effect for the string configurations obtained here. In this case, we could expect to get a 
positive intercept as given for the closed string case in [8] . We, however, do not consider 
the quantum corrections in this paper, and we concentrate on the slope parameter of 
the Regge trajectory at large E and J, where the classical approximation would be 
good. 

4 Baryon Regge trajectory and rotating vertex 

The baryon is constructed from the vertex and fundamental strings. The vertex 
is given by the D5 brane which wraps of the metric (12. 3 p (or (12. 5p ) given above 
and couples to the 5 form fluxes of the bulk through the U{1) gauge fields in it as 
shown below. Due to this dissolved fluxes, the baryon vertex could have non-trivial 
configurations. 

In order to show explicit configurations of the baryon vertex considered here, we 
parametrize the by the five angle variables, {^,^25^35^45^5}- Then the coordinate 
X^(i = 4 ~ 9) of (I22D are written as, 

= rcosO , 
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X = r sin ^ cos 62 , 

X = r sin 9 sin 6*2 cos ^3 , 

X^ = r sin 9 sin 6'2 sin ^3 cos 9^ , 

= r sin ^ sin 6*2 sin ^3 sin 9^ cos 9^ , 

X'* = r sin ^ sin 6*2 sin ^3 sin ^4 sin 6*5 , (4.1) 

and we have 77^ = I]i=4 = ''^^ sin^ 9 sin^ 6^2 . 

In general, we can consider many kinds of configurations for the D5 brane. They 
will be shown by introducing the external coordinates on the D5 brane as, r{9i) and 
X^[9j), z = 1 ~ 3, but it would be very complicated to solve the equations of motion to 
minimize the action for these general functions. Then, we restricted our analysis to a 
simple case with r{9) and (X^(6'))^ + (X^(6'))^ = p^{9). In this sense, our configurations 
adopted here are not special but simple. For more complicated configurations, they 
are remained as future works. 

The two typical configurations of baryons considered here are shown in the Figsj2] 
and E] in the three dimensional space of {X^, X^, 77}, where the distance from the origin 
is equivalent to the coordinate r, which is denoted as r = (X^)^ + (X^)^ + rf. 

As mentioned above, the D5 brane is embedded to minimize its action with the 
coordinates {r(6'), ^(6*)}, where = (X^)^ + (X^)^. Unfortunately, in the Figsj2]and 
El we can not see the coordinate p and then the extension of p{9) for the split vertex. 
The point and split vertices are discriminated here by the number of the cusp, namely 
one for the point and two for the split as shown in the Figs. [2] and El 

As for the D7 brane, its configuration is expressed by the solution given in (12. 9 p in 
the form of a tube of w =constant. We should notice here about the number of the D7 
branes. In the Fig. [21 the point vertex is shown with one cusp of D5 brane. In this case, 
we need one D7 brane at (X^, X^) = (0, w) as the end point of strings as shown in the 
figure. On the other hand, we need at least two flavor D7 branes at (X^, X^) = (0, ±w) 
for the split vertex baryon as shown in the Fig. El This implies that the number of 
the flavor quarks might be related to the configuration of the baryon vertex, especially 
to the number of the cusps on the D5 brane. This point should be examined more by 
considering more complicated configurations in the future. 

As mentioned above, both D branes are embedded as probes. Then we neglect the 
back reactions. However, we should notice that there are cusps in D5 brane due to the 
dissolved bulk flux as U{1) fields in it. On the other hand, the fundamental strings, 
which terminate on these cusps, are dissolved in the D5 brane [HI [19] as the U{1) flux. 
Then there is a smooth flow of the strings from D5 brane vertex. 

On the other hand, the embedded D7 brane has no cusp, and then strings end on 
this brane with the Diriclet conditions. This situation is the same with the case of 
mesons, which are given as strings whose two end-points are on the D7 branes with the 
same boundary condition. For the case of the baryon, however, the situation would be 
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Fig. 2: The typical configuration of point vertex(D5 brane) and flavor brane(D7 brane) in 
the space {X^, X^, r]). The point vertex has one cusp at 9 = n, and then all N^. strings 
should attach to this cusp. 




Fig. 3: The typical configurations of split vertex and two flaver branes on {X^, X^, rj). The 
split vertex has two cusp at 9 = 0, vr. uNc strings should attach to the cusp at 6* = 0, and 
Nc{l — v) strings should attach the other cusp. Two shapes of the vertex are shown. They 
depend on the boundary conditions of embedding equations, the left one with a spindle-shape 
and the right with the peanuts form in the figure. 
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slightly different since the baryon number is not zero due to the dissolved strings. As a 
result, the U{1) gauge field in the D7 brane would be a non-trivial configuration which 
reflects the non-zero charge density. This effect however could be neglected since the 
density of baryon number in the D7 brane would be negligibly small. 

In the next, we consider the vertex part, which is given by probe D5 brane. 

4.1 Rotating D5 brane vertex 

The D5-brane action is thus written by the Dirac-Born-Infeld (DBI) plus Chern-Simons 
term 

Sd5 = -T5 I d'^e-''^- det {g^b + 27ra'Fab) + T, J (2WF(2) A C(4))^ ^ , (4.2) 

where T5 = l/(^,(27r)^//) and C(4) are the brane tension and induced four form respec- 
tively. The Born-Infeld term involves the induced metric g and the U{1) worldvolume 
field strength F(2) = dA(^iy The second term is the Wess-Zumino coupling of the 
worldvolume gauge field, and it is also written as 

S = -T,J d'^ e-V-det(^ + F) + r5 J A ^(5) , 

in terms of (the pullback of) the background five-form field strength G(5) = (iC(4), 
which effectively endows the fivebrane with a U{1) charge proportional to the solid 
angle that it spans. 

Then we give the embedded configuration of the D5 brane in the given lOD back- 
ground by minimizing the action. At first, we fix its world volume as C,"" = {t, 6,62, ... , 6*5) 
For simplicity we restrict our attention to the metric fl2.12p with the ansatz, p = p{9) 
and 9 = ut. Next we consider the 5*0(5) symmetric configurations of the form p{6), 
z{9), and At{9) (with all other fields set to zero), where 9 is the polar angle in spherical 
coordinates. The action then simplifies to 




(4.3) 

where = Stt^/S is the volume of the unit four-sphere. 

The gauge field equation of motion following from this action leads 

dgD = -4sin^^, 

where the dimensionless displacement is defined as the variation of the action with 
respect to E = F^g, namely D = 5S/5Fte and S = S/T^Vt^^R^. The solution to this 
equation is 

D = D{u,9) = ^{uTx - 9) + ^ sin ^ cos ^ + sin^ 9 cos 9. (4.4) 
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Here, the integration constant v is expressed asO<z/ = fc/A^c^ Ij where k denotes 
the number of strings emerging from one of the pole of the 5*^, which is wrapped by 
the D5 brane vertex. Next, it is convenient to ehminate the gauge field in favor of D 
and Legendre transform the original Lagrangian to obtain an energy functional of the 
embedding coordinate only: 



\ 



p'J^) 1+22 




dz 



(4.5) 

V^{e)=D{y,Qf ^^m^Q , (4.6) 

where we used T-^^/^R^ = N^/ (Stt^o'), and the integral variable is changed from 6 to z. 
So 6* is a function of z in this case. Using this expression f l4.5p and (14. 6p . we can solve 
the spinning D5 brane configurations. 

In the following, we consider two cusps at 6* = and 9 = n, where the displacement 
is finite and actually we have 

D{iy,0) = ^un , D{u,n) = ^{u - l)n . (4.7) 

They should be cancelled by the fundamental strings to form a baryon. We use the 
above setting of the D5 brane for two typical baryon states, (a) One is set a.s u = 
or 1 and p = 0. We call this state as point vertex baryon. The D5 brane is seen as a 
point in our real three space since it does not extend in the direction of p. (b) Another 
configuration is obtained for < z/ < 1 and dp/dz ^ 0. In this case, the D5 brane is 
extended as a string in the direction of p, and the quarks are attached on both side of 
this extended vertex. The quarks are separated to Nci' and A''c(l — i') quarks on the 
end points of the vertex. This is called as the split vertex baryon. The configurations 
of these two types of baryonic states are shown in the Figs. [2]and[3l In the next, we 
examine the spin and the energy of these configurations. 



4.2 Spin and Energy for Baryon 
(a) Point vertex; 

In this case, the D5 brane vertex contributes only to the energy and not to the spin 
J, which is obtained from the strings attached to the D5 brane point-vertex. Then we 
have 

Jp = NJs , Ep = N,Es + Un^ip = 0) |.=o , (4.8) 

where Jg, Eg and [/ds are given by (I2.15p . fl2.16p and (14.50 . The configuration of D5 
brane for z/ = 1 is the same as the one of z/ = case, so it is sufficient to consider only 
the above case. Jp and Ep are calculated as follows. At first, solve the equation of 
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motion for z{9) by using U^^^p = 0)|i^=o- Then we find the tension of the D5 brane at 
6 = n, and this is balanced by the string part, and this is expressed by the "no-force 
condition" [23]. This condition provides the boundary condition for the equations of 
motion of strings given by ( 12.14p . It is obtained at p = and 6' = vr as, 



This condition gives the balance of the force in the z direction at the connecting point 
of the strings and the D5 brane. The condition for the force in the p direction is 
satisfied by putting the strings symmetrically around the D5 brane vertex as shown in 
the FigH 



As shown in [23] , we find two typical string configurations in the case of the point 

vertex. They are discriminated by the conditions, (i) 1/2^(0) < l/z(7r) or (ii) l/z{0) > 
l/z{n) of the D5 brane configuration. They are shown in the FigJH 

For the case of (i), the strings can not stretch so much, then this configuration 
contribute only to the small Jp and small Ep baryons. In this case, the lowest energy 
of the baryon is obtained when the strings shrink up to the D7 brane, and then the 
baryon is almost constructed by the D5 vertex. This configuration gives the lowest 
baryon mass, which is finite contrary to the lowest zero meson mass. Here we give 
the lower bound for this lowest baryon mass. It is obtained from the D5 brane which 
doesn't stretch to ^-direction {dz/d6 = 0). This is realized when the position of the D5 
brane z = is located at the same point of the D7 brane z = zdj, namely zd^ = z-qj. 
Such a solution is given [23] as follows. 





(4.9) 




Fig. 4: Rotating baryon configurations for point vertex. 





(4.10) 
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Fig. 5: The baryon trajectories with point vertex are shown for N^. = 2 5. The curves are 
given for a\.j = 0.88 (GeV-^), rUg = 0.3 (GeV), and A = R'^/a' = 10. The squares and dots 
denote the A and N baryon series respectively. 



For the above constant-point solution, we have 



«m"'^V/4j(o), jiu)= Tde^vM, (4.11^ 



where j(0) = 7.6. When we use the experimental relation, a'j^ = 0.88 (GeV~^) we find 

Ud5 = 0.S87 X NX^"^ . (4.12) 

In general, 2;d5 > ^d?, then for A > 1, we have the lowest mass bound, which exceeds 
the nucleon mass for Nc = 3. 

For large Jp and Ep, the baryons are expressed by the the configuration (ii), we can 
approximate -Ep as 

Ep ^ N,E, ~ iV,il±^Ejv/ , (4.13) 

since the vertex part Ub5{p = 0)|j/=o is negligible compared to the string part which 
is stretched long. The second equality is obtained since one of the Nc strings can be 
replaced by the half and some extra (/3) part of the meson (see (ii) of the FigJlj). Then 
we find following linear relation between Jp and Ep, 

Nc{l + 13) 2txtm 
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Then the slope of the hnear trajectory at large Ep is given as 



Regge _ Regge / . . p-n 

"iV,(l + /3)"*^ • ^^-^^^ 

This is smaller than 2/Ncy- a^^^'^, which is expected from naive quark model. Accord- 
ing to our numerical analysis, it is about 60 percent of the expected one for Nc = 3 
here. The reason of this suppression can be reduced to factor 1/(1 + /3) in equation 
fl4.14p . In the present case, we find /3 ~ 1/8 in our numerical calculation given here, 
but this factor will depends on the configurations of the rotating strings. 

On the other hand, the experimental result given by the nucleon series implies 
Qi^aae ^ g gg (^QQ^-'iy Thls can not be obtained from the above result (14.151) . Then 
we should consider another baryon configuration for higher spin states. 

(b) Split vertex; 

In the next, we consider the split vertex, and the baryon is expressed almost by the 
vertex only. In other words, the strings are absorbed into the vertex and their end 
points are attached on the two D7 branes. So the baryon configuration considered here 
is described by the D5 brane, see Fig. |H Here the two cusps of D5 brane are considered 
as the strings with zero length, then we impose the boundary condition for the cusps, 

5.p|.=.BT = 0. (4.16) 

From (14. 5p the spin is obtained as 



We remember that ( 14. 5 p originally includes the dO/dt = u, so it is further Legendre 
transformed with respect to U. Then we obtain the energy 



Ed5 = —-^^ + U = ^ I dz —e ' A^ 



duj Svr^a' J z"^ \ 1 — p^u"^ 



The energy and the spin are evaluated as in the case of the meson. Firstly, we rewrite 
the above equation (14.180 as 

Eo.=iypMS!f^, (4.19) 

MS! = :^e'"'AHz)!^^l + z^(d,ef + A-^(z){d,zr-.JvM , (4.20) 
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Fig. 6: A typical rotating baryon configuration for split vertex with minimum energy is 
show in the 3D space of (X^,X^,r). In this configuration, we can ignore F-strings which 
connect the D5 and D7 branes since their lengths are almost zero. Therefore F-strings do 
not contribute to the energy and also to the spin of the baryon. Here we notice that the two 
D7 branes at different are sitting at the same r, and the two end points of D5 brane are 
moving on the different D7 branes. 

where pi and pr denotes the end points of the D5 brane vertex. For the case of mesons, 
the end points of the string are given as pr = —pi = pb since the configuration is 
symmetric. But in the present case, the D5 brane vertex is not symmetric except for 
u = 1/2. In spite of this fact, we find a similar relation of energy and spin. 

For the case of large E05, can be approximated by the value at 6* = 6c, where 
V,y{6) is maximum. And it is given as the solution of the following equation [23] . 





nu = 9c- - sin(26'c) . 



(4.21) 




(4.22) 



on 



(4.23) 



where tm is given above Eq. (13. 6p . Then 




(4.24) 



where we used oop^ = 1 and oopL 



1. is similarly obtained as 




1 



(4.25) 



2tvtb 
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Then we could obtain the similar result to the meson trajectory with the different 
slope, 

= -i- . (4.26) 

ZTTTb 

The experimental data of a'^ lead to a similar value with that of meson, a'^ ^ ~ 1 
(1/GeV^). Then we expect 

2N 

-^sin=^(e,)~l. (4.27) 

OTT 

For Nc = 2, 3, we can not obtain a configuration for the split baryon which satisfies the 
no force condition. The reason is that one string must couple to one of the cusps of D5 
brane, and this situation is not allowed due to the no force condition [25]. Then, we 
find ^ sin3(^e) = 0.85, 1.02 for {N^, u) = (4, 1/2), (5, 2/5) respectively. In these cases, 
we can set no force condition and satisfactory slope parameter. These are assured by 
the numerical calculation given below with numerical exact solutions. 

Numerical estimation 

In order to obtain the spin and energy, we solve the equations of motion of the action 
fl4.5p . It is convenient to rewrite U in the parametrization invariant form to solve the 
equations. 



n ^ 



where dot denotes the derivative with respect to the introduced parameter s, for ex- 
ample z = dz/ds. We rewrite U as 



where 



U = ^^ldsU, (4.29) 



U = Q^z^ + zW + A^p'', Q = ^e'^/^AMl - p^cu^) V,{e) . (4.30) 



z^ 



Then the canonical momenta of qi{= {z, 6, p)) are given by pj = |^ as 



9^ A2P 



P^ = ^, Pe = z—, ,pp = A—, (4.31) 



where 



Jz^ + zW + A^p^ 
A = ^ — . (4.32) 



Then the Hamiltonian H is written as 
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H = T.P^<i^-U = Aipl + ^^+^^-Q'\ . (4.33) 
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Due to the reparametrization invariance of the system, we have the constraint H = 0. 
This constraint can be changed to a more convenient form by changing the Hamiltonian 

as 

H = 2HA, ^=^P' + f| + §-Q') ' (4-34) 
and the equations of motion are obtained for new Hamiltonian H as 




Fig. 7: The baryon trajectories with rotating D5 spHt vertex are shown for Nc = 4: (left) and 
Nc = 3 (right). The curves are given for a'^ = 0.88 (GeV-^), l/zor = 1.0 (GeV), and three 
values of A = B? /a' . The large and small dots denote the A and N baryon series respectively. 

In Fig. [71 we show our numerical results of the split baryon trajectories for the cases 
of {Nc = 4, u = 1/2} and {Nc = 3, u = 1/3}. In each case, the calculations are given 
for mq{= y/X/ {2'nzo'j)) = 0.16, 0.51, 1.04 (GeV). The slope parameter is compatible 
with the experimental spectra. However, in the case of baryon, the masses become very 
large for large A due to the vertex energy. The vertex does not disappear when the 
distance between two ends of split vertex becomes zero, and then it remains to give the 
lowest mass of the baryon. While we expect that this lowest mass would be near the 
nucleon mass, it depends on A as shown in the Fig. [3 Since the holographic approach 
is trusted for A ^ 1, it seems to be difficult to get a lowest mass comparable to the 
nucleon mass. It may be about three or four times of the nucleon mass for A ~ 5 > 1. 
We should notice that the meson trajectory is not changed by the parameter A as 
shown above. 

Another problem is that we can not get the no-force condition for the split baryon 
of A'^c = 3 [25] . The three F-strings going out from the two pole-points of split vertex 
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are separated to one and two in the case of A'^c = 3. We impose the no- force conditions 
on both sides where the F-string(s) and D5 brane couple. However, we know that on 
the side of one F-string it is impossible to balance the force due to the tension of the 
one F-string and the one of the D5 brane since the former one is always greater than 
the latter. As a result, the string shrinks up to zero length and disappears. Then this 
cusp of the D5 vertex can not be in a stable state. So we should extend to = 4: to 
have a stable split vertex. This is the reason why we show the example of numerical 
simulation only for Nc = 4: a.t the same time. 

One way to allow the baryon state of A^c = 3 and i/ = 1/3 is to introduce a pair 
of quark and anti-quark at the cusp with one F-string coupled. Then the no force 
conditions are satisfied in this case, but there appears another possible instability to 
annihilate the introduced a pair of quark and anti-quark. The energy to make this 
vertex state is however very small since the cusp is nearly at the flavor brane position, 
so we can consider this state within an energy fluctuation of quantum corrections which 
would be important at this cup point. In this sense, it is meaningful to consider the 
split vertex also for A'^c = 3 

As a result, we could say that the higher spin baryon states are realized by rotating 
vertex, which extends as a string with a definite tension. The F-strings, which are 
attached to the vertex, shrink to almost zero length and don't contribute to the spin 
of the baryon state. This picture of baryon configuration is different from the one 
considered by the naive quark model. We could find that the main ingredient of the 
baryon is not the F-string but the vertex, see Fig. [61 

5 Summary and Discussion 

Using a holographic model for mesons and baryons with higher spins, Regge behaviors 
are studied. In our model, confinement is realized due to the gauge condensate (F^) 
which determines the t ension tm of the linear potential between the quark and anti- 
quark as Tm = J{F^)/2. As a result, the slope of the meson trajectory is given by 



this quantity as a'j^ = 1/ (Try (-^^)), and then (F^) is fixed from the experimental data 
of meson spectrum. We obtain {F^) = 0.13 (GeV^) from the experimental data of p 
meson trajectory. On the other hand, the intercept of the trajectory is always zero in 
our model since the meson is approximated by a classical configuration of Nambu-Goto 
string, which is embedded in a bulk dual to the confining gauge theory. In this model, 
the string is reduced to a point when its energy approaches to zero. Then the angular 
momentum also tends to zero. This would be improved by considering the quantum 
effects, which are not considered here. 

As for the baryons, the Regge behavior is examined for two typical configurations, 
which are discriminated by the vertex. The baryon vertex is given by D5 brane which 
is introduced in the bulk background as a probe brane with U{1) flux. It wraps on 
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in the bulk. We consider two simple vertex configurations. One is given by a solution 
being extended to the fifth coordinate r as r{6). In this case, the vertex is observed as 
a point in our 4D space-time, so we call this as a point vertex. Another configuration 
is obtained as a solution of the equations of motion of two coordinates, r{9) and x{9). 
This configuration extends like a string in the plane of r — x, where x is one of our 
three space coordinates. This is called as split vertex since the vertex is observed as 
a line extended in the x direction. (See also FigsJ2] and E]) We should notice that we 
need at least two fiavor D7 branes for the split vertex baryon as shown in the Fig. |3l 
This implies that the number of the fiavor quarks might be related to the number of 
the cusps on the D5 brane baryon vertex. We should examine this point more in the 
future. 

In the case of the point vertex (see Fig. S]), spin and energy are mainly given by the 
strings. Then the slope of the baryon, a'^, is related to the meson slope as in the form 
a'^ = 2a']yj/ (iVc(l + (3)), where (3 denotes a positive factor which depends on the string 
configuration. Then, in general, a'^ < a'j^.j for > 3. As a result, we could not obtain 
the slope a^, which is comparable with the experimental data. 

On the other hand, for the split vertex case, the spin and the energy are supposed 
to be provided by the vertex only since the energy is minimized for this configuration 
(see Fig. |6]). Then we can neglect the contribution from the strings in the analysis. 
The vertex with higher spin, in this case, extends in one spatial direction x of our three 
dimensional space. Then the situation is similar to the case of mesons. That is to say 
that the slope parameter is determined by the tension of this string-like vertex when 
it extends in the x direction. This tension is in general smaller than the summation of 
the independent N^. strings, then the slope of the trajectory becomes larger than that 
of the point-vertex case. Actually we could find a slope, which is compatible with the 
experimental data, for Nc = 3. In this case, the problem is that we must adjust the 
parameters to obtain the lowest mass of the baryon compatible with the nucleon mass. 
This is performed by using unfavorably small A. Furthermore, in order to satisfy 
the no-force condition at the vertex, we should extend Nc to a large value Nc > 4. 
However, this extension moves the baryon mass to the large side. These points should 
be resolved by including quantum corrections. Another possible improvement will be 
found in other holographic model of different bulk background and probe branes. They 
are remained as future works. 
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